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Вкраплення додаткової інформації 
в частотну область цифрового сигналу 


Запропоновано метод вкраплення цифрових водяних знаків в частотно-часове представлення сигналу. 
Поворот обраних фазових складових сигналу на певний кут здійснюється за рахунок незначної зміни 
певних відліків початкового сигналу. Проведені теоретичні дослідження з метою зменшення кількості 
обчислень. Запропоновано послідовність дій для вибору оптимального значення вектору змін. 


Потреба створення універсальних стеганографічних методів захисту цифрової інфор- 
мації є однією з найбільш актуальних наукових задач сьогодення. Вкраплення цифрових 
водяних знаків (ЦВЗ) у цифрові сигнали (ЦС) є одним з методів захисту цифрової інфор- 
мації. Такий захист вважається найбільш перспективним, тому його розвиток в першу 
чергу потребує теоретичних досліджень для створення надійних методів вкраплення |. 
Оскільки факт наявності ЦВЗ не завжди приховується, головною вимогою до побудови 
стеганосистеми з ЦВЗ є стійкість до типових активних атак на стеганоконтейнер, серед 
яких згладжування, стиск, зміна формату та ін. (2). При такій постановці задачі під стій- 
кістю ЦВЗ прийнято розуміти складність видалення або псування ЦВЗ без порушення 
функціональності контейнера. Вкраплення ЦВЗ може здійснюватися як у просторовій об- 
ласті, так 1 у частотній. Частотне представлення сигналу, як правило, отримують за допо- 
могою дискретного перетворення Фур'є (ДПФ), дискретного косинус-перетворення або 
вейвлет-перетворення. 

При використанні ДПФ зазвичай інформацію вкраплюють в амплітуди сигналу |3Ї. 
За рахунок певної надмірності інформації в амплітудах таке втручання не впливає на 
функціональність сигналу. З точки зору забезпечення стійкості до активних атак краще 
вкраплювати ЦВЗ у фази сигналу, оскільки більшість методів обробки сигналу вико- 
ристовує амплітуди при обробці сигналу. Втручання у фази призводить до значного 
спотворення сигналу при оберненому перетворенні, оскільки у загальному випадку при 
зміні фази компонентів повернутися в дійсну просторову область не завжди вдається. 

У запропонованому методі ЦВЗ вкраплюється в частотно-часове представлення 
сигналу поворотом обраних фазових складових сигналу на певний кут за рахунок незнач- 
ної зміни певних відліків початкового сигналу. Функціональність сигналу при оберне- 
ному ДПФ при цьому не порушується. Для вилучення ЦВЗ до стеганоключа додаються 
номери компонентів, фази яких при вкрапленні змінювалися до потрібного значення. 

Найпростіший і водночас найдовший шлях вирішення поставленого завдання -- це 
послідовна зміна всіх відліків вхідного сигналу доти, доки фаза потрібного компонента не 
досягне шуканої величини. В даній роботі ми намагаємося зменшити кількість обчислень, 
теоретично досліджуючи, як зміни в просторовій області відображаються в області Фур'є. 

Нехай маємо дійсний дискретний сигнал є довжиною М відліків (контейнер) 


2 «|на Вуз Ї: (1) 
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Дискретне перетворення Фур'є |4| переводить його в комплексну область 
бе оз ДО 


тр Урок Я г СА росу М- І. 


Зміна відліків сигналу Є, на величини Д,, п-0,М-1 , призводить до зміни 


значень компонентів С, відповідно до виразу С, - С, ЯК, |5Ї, де 


ри з ЯЗ ян 2 


Скориставшись формулою Ейлера ехр(іф) - со5(фФ)- і51п(Ф) , довільне комплексне 


число ехр(і Ф) представимо у вигляді одиничного вектора Єр що е сов(ф)- е2 зіл(ф) 


у двовимірному ортогональному базисі (21,62), де Ф| та Є; - вектори одиничної 
довжини (рис. 1). 


е со5(Ф) Ке 


Рисунок 1 - Представлення величини ехр(ір) в ортогональному базисі (е|,Є» ) 
Введемо позначення: а, - амплітуда компонента С, (довжина вектора С,,), Фу - 
- ' 
фаза С, (кут, який утворює вектор С, з дійсною віссю); й, та (ф, - Дф,) - відпо- 


відно амплітуда та фаза С; П та 0, - відповідно амплітуда та фаза К,. З ураху- 
ванням позначень справедливі вирази 


С, з а, (б; со8ф, - ба 8іпФ, ), 


г Р Ві совіру о у б» 8іП(Фу КА 2), 
В - г, (6 со8д0, -- бо віп, ). 


Задача зводиться до пошуку вектора ВЕ, 
У векторному просторі значення доданка К, можна інтерпретувати як різницю 
векторів с, та С, , кут між якими складає ЙФ, (рис. 2). 
З трикутника, утвореного векторами С», С, та К,, за теоремою косинусів 
(п) «(а ка, ) - Зада, сов (АФ). 0) 


во... й 


і й 
-боозо- отримуємо зіп у, 2 О8іп(Л9,). 
зіПДф, 5іпу, г, 


З теореми синусів 
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З рис. 2 знаходимо 


дуаФнаФ,нУуу. (4) 
З (3) знаходиться амплітуда Гу, з (4) - фаза 0, компонента К,: 
дуяФьНаФЯ му й (5) 
Ту 
з обмеженням 

л аязп(Д л 
п «рн да АР р (6) 

2 г 2 


З урахуванням співвідношень (3), (5) 


. ш(Д 
К,ак а (з ДФ, нок зе і (7) 


Ту 


Прирівняємо праві частини (2) та (7): 


м-1 ; 
- УА, В - 3 -луехр(- фу - АФ ну,)). (8) 
пад 


Ч 
Рисунок 2 - Вектор різниці К, -С,- С, 

Рівняння (3) обмежує амплітуду п,. Для того щоб зміни початкових відліків 
сигналу не були помітні в області Фур'є, а, та а, не повинні сильно відрізнятися. 
Позначимо (а, - ,) - д , тоді з (3) 

2 
2 2 б 5 
(5 - (а,) Древо 2 (9) 
а, а, 

З (9) видно, що для заданих й, та ДФ, обмеження по б приводить до обмежен- 

ня по Р,. Оскільки за умовою задачі С) відомі, відповідно відомі й Ф, та відомі 


" 
бажані Ф,, тому ЙФ, теж відомі. За допомогою (5) знаходимо 0, при відомих /у. 
Потрібно мінімізувати зміни амплітуд компонентів сигналу. Крім того, рівнян- 


ня (9) обмежує зміну кута вектора С). За умови, що довжина вектора К, менша довжи- 
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ни вектора С,, можна оцінити найбільшу величину кута ЛФ, (рис. 3). Якщо навколо 
кінця вектора С, окреслити коло радіусом 7, та з початку вектора С, провести дотич- 
ну до цього кола, тоді кут між дотичною та самим вектором С, буде максимально 
Ж М 


можливим кутом ЛоФ,,. Враховуючи, що 5іп АФ - т--1,, отримуємо ЙФутах - аГС8іП 7 


У У 


Бажано, щоб У, - 7/2. Тоді з (4) отримуємо 0, - Ф, - ДФ, - 7/2. На практиці 
ми не завжди можемо вносити зміни у відліки початкового сигналу так, щоб досягти 
такого значення У; , але завжди треба пам'ятати, що саме такий кут є бажаним для 
зменшення змін відліків вхідного сигналу. 


2, 


8, 
с і. 


| 


Рисунок 3 - Ілюстрація максимального значення кута ДФ, 


Послідовність дій для вибору оптимального значення К, може бути наступною. 
При відомих у та м визначити амплітуду та фазу у-го компонента, вирахувати ДФ, . 


Відповідно до рис. 3 від початку вектора С, під кутом ЛФ, провести промінь. З кінця 
вектора С, опустити на цей промінь перпендикуляр. Напрямок даного перпендикуляра 


буде найбільш оптимальним напрямком вектора змін К,. У випадку, коли змінюється 
значення тільки одного відліку початкового сигналу, потрібно змінювати відлік з тим 


Жпіпу ран - 
було найбільш наближеним до оптимального напрямку К, 


номером п , щоб - 


В роботі (51 розглянуто випадок, коли амплітуди С, та С у збігаються. Для нього 
була доведена наступна теорема. 
Теорема. Для зміни фази С, - аехр(їФ) на величину Д и необхідно додати до 


С, величину а І (о з дея гетяі |. : 


а(ехр(іф А, )) - а(ехр(іф)) - Завіп пе ак (о з зеія 5)| | ( 10) 


Наслідок. З теореми випливає, що змінити фазу Ж, - аехр(іф) на величину ДЛ, 
можна, вносячи такі зміни в просторову область, щоб 


зазін о (оч зезяї аа З али (11) 
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Повернемося до загального випадку, коли змінюються і фаза 1 амплітуда у -го 
компонента в області Фур'є. Розглянемо найпростіші випадки, коли змінюється один 
або два відліки вхідного сигналу. Нашим завданням є зміна Ф, на /, -Ф, - АФ, 


(далі м, - у ), з найменшими змінами відліків вхідного сигналу, тобто з найменши- 
ми значеннями А,, п-0,М-1 . 


Випадок 1. Зміни вносяться лише до одного відліку сигналу є в (1). Нехай 


араяро А (А 2), тоді 


д тку 
Отже, п - ГІ 9,2- х з обмеженням (6). 
Потрібно визначити К та У такі, щоб 4 було найменшим за абсолютною 
величиною. З (5) отримуємо 


. | Ма, зіп (АФ, ) 2лку 
агсвіп г. 
Д М 


Ма,,він(АФу,) |. пк 
, або У УЛ іп и 
и - -2Еачу 


або розглядаючи Д як функцію від К 
Ма, зі ( АФ, ) 


ПА ПСИ ОВО 
"Бо. у" з 
зщ| - ооо 

М 


(Ю) може бути знайдена розв'язанням наступного 


ЛЮ з 


Шукана величина Д 


тіп 
рівняння 
тіп Ма, зі (АФ, ) 


Д ії (Ю) , при на оби (12) 
нео - -м) М 


З (12) видно, що потрібним є значення /, для якого зі ((-27ку / М) -м)| досягає 


найбільшого значення за абсолютною величиною, тобто кут (тку /Муч у) наближе- 


ний до кута їл/2. 
Розглянемо найпростіший приклад: нехай сигнал 2 -|2є0, 71,52 у-1| Є ЛінійнимМ 


2; З ад З ібо , довжина блока М «8. Як показали дослідження у роботі |5|, найменше зна- 
чення амплітуди в області Фур'є буде мати величина Сз, тому обираємо у - 3. На рис. 4 
показані напрямки всіх складових вектора 5, Також на рис. 4 показано, як формується 
вектор С,, його довжина та напрямок. З цього ж рисунка видно, що для невеликих 
значень Дф, найкращим рішенням буде зміна п'ятого відліку (для від'ємного ДФ, ) або 
шостого (для додатного Д9, ). Рівняння (11) дає той самий результат. 

Випадок 2. Зміни вносяться до двох відліків сигналу є в (1). Нехай Р зяїча, 


" 


Єт 7 Єт Я 2, причому КП, В області Фур'є ці зміни відіб'ються на значеннях 
З ен | 
з ехрі - і 
М 
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Скориставшись формулою Ейлера та тригонометричними рівняннями, можна 
виразити В, одним комплексним числом 


А. ліку 2лту (зидле -  дяту 
К, ---| со з со5 й 8іп віп - 
М М М М М 


а Р АН п т) зт пров ст М 
М М М М М 


з ек і б пі во з т)- тей з п?) че 
М М М М 


2А 2лу 2лу 
2---с085---(К -т)ехрі - і-- (КЕ . 
З бу (к-т) | б ( п) 


Отже, В, - п, ехр(-і0,), де п, - зовн -т), в,з А -т) з обмеженням (б). 


ь - 3 
е1 


з 3 
т Єб 
Рисунок 4 - Формування вектора Сз для випадку 2 
Проводячи міркування, подібні до розглянутих вище для першого випадку, роз- 
ри . | Длм 
глядаючи ДА як функцію від К , отримуємо при "Ба (к-т)-фу - що -0 


тіп Ма, іп (АФ, ) 


ДЛ (К)с й 
Я | ) кед МІ : лу 
2 зп р УЮ Дрова 


(13) 
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У випадку 2, - 8, ЗД, є, 7 Є, - А, провівши перетворення, аналогічні по- 
переднім, отримуємо 


р і ча Р сло) є ен) 


Проведені дослідження показали, що в загальному випадку, якщо всі відліки 
сигналу змінюються на величину ДА (або лишаються без змін), амплітуда К, про- 


Жпіпу 
М 


оптимальне значення ЛА з виразу (5), перебираючи тільки номери відліків, які зазнають 
змін під час вкраплення. За рахунок цього обчислювальна складність пошуку розв'язку 
поставленої задачі значно зменшується. В частинних випадках, розглянутих в роботі, 
вдалося отримати аналітичні рівняння для знаходження оптимального розв'язку постав- 
леної задачі. 

Зауважимо, що, з точки зору забезпечення функціональності сигналу, необхідно 
забезпечити відсутність сплеску амплітуд, особливо для низьких частот. Рівняння (9) 
дає зв'язок між величинами Д та 5, що дозволяє контролювати сплески амплітуд. 


М-1 
. А А З 
порційна м тобто К, од зб . Це означає, що ми зможемо знайти 
под 
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Л.Л. Никитенко, О.Ю. Никитина 

Встраийвание дополнительной информации в частотную область цифрового сигнала 

Предложен метод встраивания цифровьгх водяньх знаков в частотно-временное представлениє сигнала. 
Поворот вьюбранньх фазовьїх составляющих на определенньй угол осуществляєтся за счет незначительного 
изменения некоторьїх отечетов входного сигнала. Проведень теоретическиє исследования с целью 
уменьшения обьема вьтчислений. Предложена последовательность действий для вьтбора оптимального 
значения вектора изменений. 


І. МіКйепко, О. Мікйийїпа 

Тие 5іде Піїогтабйоп Нідіпе Піо (пе Ріосна! 5іспа! Егедпепсу Ротаїп 

Ніфпє тефод ої (Пе аїеїаї ууаїегплагк їпіо тпе-Їтедиепсу 5іспа! поіайоп м/а5 ргорозед. ТБе 5еіссівад 
рбазе їегта гобайоп оп Пе ргебеї апеїе аз саггтуїпо, оці мїб іпзідпійсапі спапее5 ої Не 5оте ої Ше іприї 
віспа! соптропепіз. Тре Шеогейїса! ге5еагсі е5 ууеге пладе улів Ше єоа! їо гедисе согарийтеє. ТРе ууогкПом/ 
умав ргорозеа іо аке сбапєе уесіог Без! уаїце. 


Стаття надійшла до редакції 06.07.2010. 
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